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ABSTRACT 

Several recent studies find that 10 - 50% of morphologically selected field early-type galaxies at redshifts z <, 1 
have blue colors indicative of recent star formation. Such "blue spheroids" might be massive early-type galaxies 
with active star formation, perhaps induced by recent merger events. Alternatively, they could be starbursting, 
low-mass spheroids. To distinguish between these two choices, we have selected 10 "Blue Spheroid Candidates" 
f\ (hereafter, BSCs) from a quantitatively selected E/SO sample to study their properties, including kinematics from 
Keck spectra obtained as part of the DEEP Groth Strip Survey (GSS). Most BSCs (70%) turn out to belong to two 
broad categories, while the remaining objects are likely to be misclassified objects. Type-1 BSCs have underlying 
red stellar components with bluer inner components. Type-2 BSCs do not show an obvious sign of the underlying 
red stellar component, and their overall colors are quite blue ((U-B) rest < 0). Both Type-1 and Type-2 BSCs have 
internal velocity dispersions measured from emission lines a < 80 km sec" 1 and estimated dynamical masses of 
only a few x 10 10 Mq or less. For Type-1 BSCs, we estimate a of the red component using the fundamental plane 
relation of distant field absorption-line galaxies and find that these a estimates are similar to the a measured from 
emission lines. Overall, we conclude that our Type-1 and Type-2 BSCs are more likely to be star-forming low 
mass spheroids than star-forming, massive, early-type galaxies. 

8 Throughout this paper, "spheroids" are defined as featureless galaxies which include Es, SOs, and dEs, while bulges of spiral galaxies are 
not included. 

Subject headings: cosmology:observations — galaxies:dwarf — galaxies: elliptical and lenticular, cD — 

galaxies: evolution — galaxies: formation — galaxies: kinematics and dynamics — galaxies: 
high-redshift 



1. INTRODUCTION 

When two gaseous galaxies merge, the gas sinks towards the 
center of the merger product, creating an intense burst of new 
stars, and the merger product has the surface brightness pro- 
file dominated by the r 1 / 4 law (e.g., Mihos & Hernquist 1994). 
Thus, the detection of many smooth, highly-concentrated, blue 
star-forming galaxies would provide strong support for the for- 
mation of early-type galaxies via mergers at z < 1, as predicted 
by some hierarchical merger models (e.g., Baugh et al. 1996; 
Kauffmann et al. 1996). 

Recently, several groups have identified blue galaxies at z <, 1 
that appear to be early-type galaxies such as E or SOs (Im et al. 
2001; Menanteau et al. 1999; Schade et al. 1999; Abraham et 
al. 1999; Franceschini et al. 1998). Their surface brightness 
profiles show a significant r 1 / 4 component typical of nearby E 
and SO galaxies, typically with bulge-to-total light ratios (B/T) 
greater than 0.4 (Im et al. 2001). Some of these galaxies are 
found to have [Oil] emission lines, indicative of recent star for- 
mation (Koo et al. 1996; Schade et al. 1999). Also, others 
find younger, blue stellar components on top of underlying old, 
red stellar populations (Menanteau et al. 2000; Abraham et 
al. 1999). These morphological features, as well as their blue 
colors and the rather high luminosities of at least some of the 
blue spheroidal galaxies, suggest that they may be direct pro- 



genitors of the present-day typical, massive early-type galaxies 
(> 10 1 'Mq or a > 150 km/sec) which are undergoing strong 
star formation after merging. 

The blue spheroidal galaxies, however, could instead be 
low-mass, starbursting dwarf galaxies such as dwarf ellip- 
ticals/spheroids (hereafter, dE following Mihalas & Binney 
1981; e.g., Im et al. 1995a, 1995b; Babul & Ferguson 1996; 
Driver et al. 1996). The strong star formation activity at the 
center might create excess light over the exponential surface 
brightness profile commonly found in dwarf spheroids (Lin & 
Faber 1983; Gallagher & Wyse 1994), which could be inter- 
preted as a significant bulge component. In fact, some blue 
spheroidal galaxies resemble distant compact galaxies (CGs; 
Koo et al. 1994, 1996; Guzman et al. 1996, 1997, 1998; 
Phillips et al. 1997). The previous works on CGs suggest 
that at least some CGs are likely to be progenitors of low-mass 
spheroids (<; 1O 1O M ) today. 

As a part of the DEEP Groth Strip Survey (GSS; Vogt et al. 
2001a; Koo et al. 1996), we have identified 10 BSCs from HST 
images, and we have measured their emission line widths to 
obtain mass estimates. In this paper, we present evidence that, 
based on their kinematic and structural properties, the majority 
of BSCs are likely to be low-mass spheroidal galaxies undergo- 
ing strong starbursts, rather than massive star-forming E/SOs. 
Throughout the paper, the rest-frame quantities and model pre- 
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dictions are calculated assuming fl m = 0.2, A = and h = 0.7 
for the cosmological parameters, where the Hubble constant is 
normalized as h = Hq / 100 km/sec. 

2. DATA 
2.1. Sample 

The BSCs are drawn from a set of 28 contiguous HST 
WFPC2 fields, the so called "Groth Strip" (Groth et al. 1994; 
Rhodes et al. 2000). A significant portion of the Groth strip 
was spectroscopically observed with the Keck 10-m telescope 
by the Deep Extragalactic Evolutionary Probe (DEEP; Koo et 
al. 1996; Vogt et al. 2001a; Simard et al. 2001; Phillips et al. 
2001). 

A combination of HST images and Keck spectra enables us 
to select BSCs as described below. We emphasize that our se- 
lection criteria for BSCs are quite different from those of the 
CGs in Phillips et al. (1997) and Guzman et al. (1997), and that 
our goal is to select galaxies which morphologically resemble 
local Es and SOs but with blue colors. 

The first step in our procedure is to fit each HST galaxy im- 
age using the 2-dimensional bulge+disk surface brightness fit- 
ting algorithm GIM2D (Simard et al. 1999, 2001; Marleau & 
Simard 1998). The main purpose of fitting HST galaxy images 
with GIM2D is to derive structural and morphological parame- 
ters that can be used for quantitative selection of E/S0 galaxies 
(Im et al. 2001). The surface brightness profile of the bulge 
component is fitted with either a de Vaucouleurs (or r 1 ' 4 ) profile 
or an exponential profile. The de Vaucouleurs profile is known 
to fit the surface brightness profile of ellipticals and bulges of 
luminous spiral galaxies (de Vaucouleurs 1948). Bulges of less 
luminous, late-type spirals are known to be fitted well with an 
exponential profile (e.g., Andredakis & Sanders 1994; de Jong 
1996). For the component fitted as a bulge, we will use the 
term "photo-bulge", since this component is mathematically de- 
fined by the surface photometry, though in some cases may not 
correspond to physically-defined "bulges" of nearby galaxies 
in terms of stellar population and kinematics. For the sample 
selection of BSCs, we will only use a de Vaucouleurs profile 
for the photo-bulge component. The surface brightness pro- 
file of the disk component is fitted with an exponential profile, 
which is known to fit the SB of disks of spiral galaxies. In anal- 
ogy to the photo-bulges, we will use the term "photo-disk" for 
this component. GIM2D measures model-fit structural parame- 
ters which include the total magnitude; bulge-to-total light ratio 
(B/T); disk scale length (r^); bulge effective radius (r e ); galaxy 
half-light radius (r/,/); position angles and ellipticities for both 
the bulge and disk components, and errors associated with each 
quantity. Centers of both the photo-bulge and photo-disk com- 
ponents are assumed to be the same, and fitted along with the 
other structural parameters. This procedure of surface bright- 
ness fitting is similar to that of Ratnatunga, Griffiths, and Os- 
trander (1999), and Schade et al. (1995). GIM2D searches for a 
parameter set which gives the maximum likelihood value using 
Monte-Carlo sampling of the parameter space. Then, GIM2D 
Monte-Carlo samples the region around the maximum likeli- 
hood, and the median value of each of the parameters is adopted 
as the model-fit parameter; the errors are estimated at the 68% 
confidence level of the statistical distribution of the sampled 
parameter sets (Simard et al. 2001). The errors may be un- 
derestimated for two reasons: (1) PSF modeling errors are not 
included; (2) the sky levels are not fitted along with other pa- 
rameters, thus the confidence intervals do not include errors on 



the sky level estimate. Our tests with simulations show, how- 
ever, that these effects do not severely underestimate errors (see 
Simard et al. 2001). The assumed luminosity profiles can also 
affect the output parameters (see section 3.2). In general, V and 
/ band images are fitted with GIM2D separately, which we call 
as "separate fit", as opposed to "simultaneous fit" where V and 
/ band images are fitted simultaneously to provide better con- 
straints on the model-fit parameters such as V - 1 (see section 
3.2). 

In the next step, early-type galaxy candidates are selected 
based on two quantitative morphological parameters, the bulge- 
to-total light ratio (B/T) and the residual parameter (R), which 
are derived from GIM2D using the "separate" fits. We use 
the "separate" fits for the sample selection since our early-type 
galaxy candidates come from the early-type galaxy sample of 
Im et al. (2001), which is based on the "separate" fit. The quan- 
tity B/T measures the ratio of the luminosity in the photo-bulge 
component to the total luminosity contained in both photo- 
bulge and photo-disk components. Thus, the B/T describes 
how prominent the photo-bulge component is, and typical val- 
ues are B/T > 0.5 and B/T > 0.3 for visually classified local 
Es and SOs, respectively (Im et al. 2001; Kent 1985; Scorza et 
al. 1998). Note that the B/T cut alone is not very effective at 
isolating E/SOs. Some late-type galaxies have B/T > 0.3 - 0.4, 
and about 30% of objects selected with the B/T > 0.3 - 0.4 cut 
are indeed non-E/SO galaxies (Im et al. 2001; Kent 1985). For 
that reason, we use another parameter - the residual parameter 
(R) - to exclude non-E/SO galaxies with large B/T values. 

The residual parameter R describes how symmetric and 
prominent the morphological substructures are (such as spiral 
arms and HII regions). Smooth-looking, symmetric galaxies 
like E/SOs generally have R < 0.08. For a detailed descrip- 
tion of the R parameter, see Im et al. (2001) and Schade et 
al. (1995). Im et al. (2001) select "Quantitatively-Selected 
E/SOs" (hereafter, QS-E/SOs) with the criteria B/T > 0.4 and 
R <J 0.05 - 0.08, depending on the size and apparent magnitude 
of the object. Tests on simulated local galaxy images show that 
these criteria select E/SOs well (Im et al. 2001). The R cut 
varies as follows: R < 0.08 when I < 21.0, R < 0.06 when 
21.0 <I < 21.5, and R < 0.05 when 21.5 <I < 22. (/-band 
magnitude is defined as the total model magnitude derived from 
the surface brightness fit; note that this total magnitude will be 
used throughout the paper, unless noted otherwise). The R cuts 
are chosen to minimize the contamination of the QS-E/S0 sam- 
ple by late-type galaxies when the S/N is low (S/N < 50). We 
use the same criteria to choose QS-E/SOs in this paper, which 
are identical to QS-E/SOs in Im et al. (2001). 

As a next step, the spectroscopic redshift-color diagram is 
used to identify unusually blue QS-E/SOs as BSCs, where the 
color (V-/) is also derived from the model total magnitudes. 
Fig. |] shows the (V-I) vs. spectroscopic redshift diagram 
of a magnitude limited sample (/ < 22) of 262 galaxies in the 
Groth strip. Note that the identification of spectroscopic red- 
shifts is virtually complete to this magnitude limit for the GSS 
(Phillips et al. 2001). Also shown in Fig. [I] are lines of plausi- 
ble color ranges for passively evolving, old stellar populations 
computed with the 1996 version of the spectral synthesis model 
of Bruzual & Chariot (1993; hereafter, GISSEL96) - the upper 
dashed line is created by assuming a 0. 1 Gyr burst with 250% 
solar metallicity, Salpeter IMF, and formation redshift z j or = 11. 
The lower line assumes the same model with lower metallicity 
(40% solar). To allow for magnitude errors, we add/subtract 
0.15 magnitude in (V-I) color to the upper/lower lines. Im et 
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al. (2001) find that most QS-E/SOs (marked with squares) form 
a red envelope in the (V-I) vs. redshift diagram. However, 
six out of 44 QS-E/SOs (marked with thick triangles) have col- 
ors that are bluer than the expected color range of the passive 
evolution model; we call these QS-E/SOs "good" BSCs. 

In addition to these "good" BSCs, we have examined im- 
ages of galaxies that would have been selected as BSCs if the 
morphological selection criteria were slightly less restrictive 
(B/T > 0.35 or R <j 0.09). This visual inspection reveals 4 
additional candidates, and we call them "possible" BSCs (thin 
triangles in Fig. [I]). Fig. || shows /-band images of the 6 "good" 
BSCs and 4 "possible" BSCs. By constructing an azimuthally- 
averaged surface-brightness (SB) profile using ellipsoidal aper- 
tures and plotting it along the major axis, we confirm that 
our BSCs show surface brightness profiles which differ signif- 
icantly from a simple exponential law, but that are well fitted 
by the addition of a significant photo-bulge component. This 
and other important points regarding the colors and structures 
of the underlying stellar populations are discussed in more de- 
tail in Section 3.2. 

Fig. H shows the spectroscopic redshifts vs. / for "good" 
BSCs (thick triangles) and "possible" BSCs (thin triangles), 
compared with red QS-E/SOs (squares) and all other types of 
galaxies (dots) in the GSS. Below z = 0.5, there tends to be a di- 
chotomy in the sense that red QS-E/SOs are intrinsically bright 
while BSCs are faint. Their faint apparent magnitudes suggest 
that low-z BSCs are unlikely to evolve into massive E/SOs to- 
day. In contrast, above z = 0.5, BSCs approach L* in luminos- 
ity. However, the apparent brightness of z > 0.5 BSCs does 
not necessarily mean that they will evolve into the present-day 
L* E/SOs, as the luminosity dimming for these galaxies may be 
more significant than that of the old stellar populations which 
are assumed in lines of Fig.^]- thus, these L* high-z BSCs 
may fade significantly into low-luminosity systems today. 

To improve constraints on the local counterparts of BSCs, we 
need to study quantities which are independent of the luminos- 
ity evolution and are more direct tracers of the mass, e.g., the 
internal velocity dispersion a. The next section describes our 
measurements of the internal velocity dispersion, a, for these 
BSCs. 

2.2. Spectroscopic observations and line width measurements 

Spectra of BSCs were obtained as a part of the DEEP sur- 
vey during the 1996-1999 Keck observing runs (Phillips et al. 
2001) using the Low Resolution Imaging Spectrograph (LRIS; 
Oke et al. 1995). A 900 lines/mm grating was used to cover 
the blue part of the spectrum (~ 5000 - 6500 A) and a 600 
lines/mm grating was used for the red part (~ 6500 - 8500 A); 
typical exposures were ~ 1 hour with each grating. The spec- 
tral resolution (FWHM^ 2 A and 4 A for the blue and red 
parts, respectively, or Act ~ 80 km sec -1 with Act defined as a 
of the instrumental resolution) allows us to measure a as low 
as ~ 30 km sec" 1 . For one BSC (092_1339), we also obtained 
a 15-min high resolution spectrum (Act ~ 17 km sec" 1 ) using 
the Echellette Spectroscopic Imager (ESI; Sheinis, et al. 2000) 
during May, 2000, at Keck. 

To measure the line widths, we fit Gaussian profiles to strong 
S/N emission lines following the method described in Phillips 
et al. (1997) and Guzman et al. (1997). For the [Oil] doublet 
at 3727 A, we used a double Gaussian profile with a fixed AA 
between both peaks. Each velocity width was corrected for the 
instrumental resolution. Finally, we take the variance-weighted 
mean value of the velocity widths measured from different lines 



as our adopted a. 

Estimating true errors of velocity widths is not a trivial task, 
as uncertainties in velocity widths include various factors be- 
yond the formal measurement error. This is particularly true 
in estimating the instrumental correction, which varies with 
galaxy extent, seeing disk size, and misalignment of the target 
with the slit. In addition, the non-uniform distribution of the 
emission and non-Gaussian velocity profiles contribute to the 
intrinsic error. Simulations with realistic parameters indicate a 
"robust" velocity error of about one-third the instrumental line 
width, a value which is supported by repeat observations of the 
same objects, and by line width measurements using different 
algorithms (Phillips et al. 2001; Guzman et al. 1997). For 
LRIS, this corresponds to ~ 30 km sec" 1 depending on the line 
and redshift of the object. Hence, we adopt l-er error, 8a = 30 
km sec" 1 , as a good (conservative) error estimate (see Phillips 
et al. 2001 for a full discussion). For the ESI data, we adopt 
8a = 6 km sec" 1 , which is simply one-third of the ESI resolu- 
tion. The a values can be found in Table 1 . 

In Fig. [|, we show spectra centered on the [Oil] emission 
line for 3 representative BSCs. Interestingly, an inspection of 
the 2-D spectra of two "good" BSCs (092_4957 and 294_2078) 
show tilted emission lines, indicating that they are rotating ob- 
jects (e.g, disks). One "good" BSC (212_1030) has a weak 
[Oil] emission line, but it is not clear whether the emission line 
comes from the object itself or from its blue neighbors, at least 
one of which is at the same redshift. Finally, one "possible" 
BSC (082_5252) has a redshift taken from the literature (Lilly 
et al. 1995), so we do not have a measured a for this object. 

3. RESULTS 

In this section, we present our two main findings about BSCs. 
In section 3. 1, we will show that the dynamical masses of BSCs 
are similar to those of nearby low-mass spheroidal galaxies, as- 
suming that the dynamical mass estimated from the emission 
line widths is not severely underestimated (less than a factor of 
~ 2). Then, in section 3.2, we show that 5-6 out of 10 BSCs 
seem to have an underlying old stellar population, and that their 
blue overall colors originate from localized regions containing 
young stars. In the same section, we also discuss a plausible 
shape of the surface brightness profile as a means to explore the 
structure of the underlying stellar population. 

3.1. Mass of BSCs 

In Table 1, we present the basic data for BSCs. Note that 
structural parameter values in Table 1 come from the "separate" 
fit described in Section 2. We estimate the dynamical mass us- 
ing the virial theorem for spherical r 1 / 4 law systems (Poveda 
1958; Young 1976; Illingworth 1976): 

M dyn ~ 3al D r 1/2 /G, 

where a^o is the central velocity dispersion of the system in real 
3-D space, r\/2 is the projected half-mass radius, and G is the 
gravitational constant. We assume that stellar mass traces light, 
i.e., r 1/ / 2 — ri,i.cir, where r;,/ cir is the projected circular aperture 
half-light radius of the stellar system in Table 1 (roughly equal 
to the major-axis half light radius, if the ellipticity, e, is small). 
Note that we use r/,; cl >, instead of ri,i. ma j which is defined along 
the major axis and given as a default output from GIM2D, in 
order to facilitate the comparison with local samples for which 
rhi,maj is not readily available. The line of sight velocity disper- 
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sion a would be cr^o / V3 in an isotropic velocity field. There- 
fore, the total dynamical mass of the system is roughly, 

M dyn ~9a 2 r hLcir /G. (1) 

This formula is similar to that adopted by Phillips et al. 
(1997). 

Fig. ||, shows measured a vs. r;,; ci > for all BSCs (triangles). 
Also plotted in Fig. |B| are regions occupied by different galaxy 
types at z = 0, taken from Phillips et al. (1997) and Guzman et 
al. (1998). Dashed lines are drawn to indicate constant mass. In 
the same plot, squares indicate r/,/ „> vs. a of distant, red, QS- 
E/SOs (/ < 22) with spectroscopic redshifts in Im et al. (2001). 
For the Im et al. (2001) sample, we use measured <j a bs when 
they are available (Gebhardt et al. 2001), and <jfp (see Section 
4.1) otherwise. Fig. |5] shows that no BSCs resemble present- 
day, massive, early-type galaxies or distant red, QS-E/SOs. The 
majority of BSCs (5 out of 9) lie within or on the border of 
the polygon for local dEs, suggesting a possible link. BSCs 
have a < 80 km/sec and M < 10 10 Mq . Typical, massive, lumi- 
nous E/SOs (L ;> L*) have a £ 150 km/sec and M £ 10 n M Q . 
Thus, BSCs are not likely to evolve into massive E/SOs today 
under the assumption that the measured a reflect the dynamical 
galaxy mass. Studies exist which support this assumption (e.g., 
Telles & Terlevich 1993; Kobulnicky & Gebhardt 2000 and ref- 
erences therein), but future studies with more extended samples 
may be needed to establish the link between a measured from 
emission lines (hereafter, a em ) and dynamical mass of galax- 
ies - these points are discussed in section 4.1, where we argue 
that our conclusion here is not likely to change even if our a em 
values somewhat underestimate Mj v „. 

Note that CGs studied by Phillips et al. (1997) and Guzman 
et al. (1996) also lie in a similar region where BSCs are (see 
Fig.9 of Phillips et al. 1997). From the B/T and R values 
which we derived separately for CGs, we find that roughly 10% 
of CGs at / < 22 in the Phillips et al. (1997) and Guzman 
et al. (1996) sample can be classified as BSCs. The number 
of BSCs found in the CG sample per WFPC2 frame (0.25 per 
frame) matches well with the number of BSCs found here. Con- 
versely, we can apply the selection criteria of CGs in Phillips et 
al. (1997) to our sample, and we find that 8 out of 10 BSCs are 
CGs. Thus, we conclude that most BSCs belong to a subset of 
CGs that are smooth and photo-bulge-dominated. 

A more detailed description of each BSC in Fig.^J is pre- 
sented below. One BSC (062_6465; #7) lies in the region of 
Irr galaxies, but not too far from the locus of dEs. Two BSCs 
(294_2078, 212_1030; #6 and #5 respectively) are found in the 
spiral galaxy polygon. As noted above, the spectrum of one 
of them (294_2078) shows tilted emission lines indicative of a 
rotating system. Furthermore, the measured B/T value has a 
relatively large error (e.g., B/T = 0.59±[] [4 from Table 1), and 
this galaxy is a highly-inclined system. Thus, BSC 294_2078 
is probably an early-type spiral with a significant disk compo- 
nent or a later-type galaxy. BSC 212_1030 is an absorption-line 
object with a possible weak [Oil] emission line. The blue com- 
panion on the left in Fig. 2 has a redshift virtually identical to 
212_1030 (z = 0.88), and is possibly physically associated. The 
projected separation between both objects is ~ 1.2 arcsec (6 /T 1 
kpc). The blue galaxies on the right are foreground objects with 
redshifts z — 0.24, much lower than the redshift of 212_1030, 

6 IRAF is distributed by the National Optical Astronomy Observatories, under 

7 Object pixels are defined as a contiguous region of pixels that have values > 
1.5 pixels; see Simard et al. 2001 for more details. 



and so are not physically associated. Unfortunately, the spec- 
trum of 212_1030 is too noisy to measure the absorption-line 
velocity width. Its nuclear color measured with an aperture size 
of 0.5 arcsec is V -I = 1.9, much redder than the total color 
V -I = 1 .58 from the model fit. We also note that model-fit col- 
ors of this object using different procedures are V -I ~ 1.7 - 
1.8. Therefore, we suspect that the model-fit color of 212_1030 
is somehow contaminated by its blue neighbors, and that this 
object is not a genuine BSC, but rather is more likely to be a 
normal, red, early-type galaxy. 

Another BSC (092_4957; #2), which lies close to the bor- 
der of dEs and spirals, has a rotation curve with an inclination- 
corrected rotational velocity, V max ~ 130 km/sec (Vogt 2000, 
private communication). The velocity width of a ~ 60 km 
sec -1 , is in agreement with the V max value considering the con- 
version factor of l f l . in V max ~ l f l , x atom, which typically 
amounts to ~ 2 (Kobulnicky & Gebhardt 2000). The disk com- 
ponent of this object is found to lie on the local TF relation, a 
strong indication that this object may be a relatively low-mass 
spiral galaxy (Vogt et al. 2001b). 

3.2. Underlying stellar population 

The accretion of gas-rich satellite galaxies or gas infall to- 
ward the centers of E/S0 galaxies may provide fuel which will 
trigger intense star formation. Young stars in the satellite galax- 
ies can also fall into the centers. The existence of relatively- 
young stellar populations among local E/S0 galaxies suggests 
that such a sprinkling of young stars over an old, underly- 
ing, stellar population has occurred (e.g., Chariot & Silk 1994; 
Trager et al. 2000). If BSCs are old galaxies sprinkled with 
young stars near their center, we should find that the outskirts 
are red, while the inner parts are blue. Fig. ^ shows the rest- 
frame U-B color gradient of BSCs along their major axis out 
to ~ 3 ryujnaj, as well as a 2-dimensional color map. To con- 
struct the 2-D color gradient plot, we used the IRAF 6 task, 
"ellipse" with a fixed object center, an ellipticity and a position 
angle which all correspond to the values derived for the photo- 
bulge component of the "separate" 2-D fit. Background sky 
levels were estimated using a region which is at least 5 pixels 
away from the object pixels 1 . In the process, observed V -I 
values are converted to the rest-frame U-B using the empirical 
^-correction which is presented in Gebhardt et al. (2001). 

To construct the 2-D color map, we used the following pro- 
cedure: First, both the V- and /-band images of each object 
are background subtracted, and smoothed with Gaussian filter 
with (7=1. This smoothing of the image with the Gaussian fil- 
ter is necessary to enhance the global color gradients. In the 
next step, each pixel value is converted into magnitude units. 
Then, the /-band image is subtracted from the V-band image to 
produce a color image. Finally, the V — I values are converted 
to the rest-frame U-B color. The expected U-B values for a 
passively evolving stellar population after a 0. 1 Gyr burst are 
U-B = 0.1 - 0.2 at the age of 1 Gyr and U-B = 0.2 - 0.5 
at the age of 5 Gyr, for the assumed metallicities of 40% and 
100% solar values, respectively (GISSEL96). 

Fig. ^ suggests that there are three categories of BSCs, as 
described below. 

Type-1: Five, possibly six, out of the 10 BSCs show a 
tendency to have blue inner regions and red outer regions, 

□ntract to A.U.R.A. Inc. 

.5<TM g after convolving the image using a 3 X 3 Gaussian kernel with a FWHM of 



Blue Spheroidal Galaxies 



5 



suggesting that these BSCs have underlying old stellar popu- 
lations (092_1339, 103_6061, 172_5049,062_6465,082_5252, 
and possibly 153_6056). Red outer regions are as red as the ex- 
pected colors of passively-evolving stellar populations (U-B a 

0. 1 ~ 0.5), and redder than colors of local Sbc's (U-B~ 0.0). 
Thus, these BSCs might be old galaxies with star formation oc- 
curring more actively near their centers; similar objects have 
been discussed in Guzman et al. (1998), Abraham et al. (1999) 
and Menanteau et al. (2000). Note that Jansen et al. (2000) find 
that 50% of the local, low-luminosity galaxies have the bluest 
colors at their centers - Type-1 BSCs may be comparable to 
such local low-luminosity galaxies. 

Type-2: One, or possibly two BSCs (273_4427 and possibly 
153_6056) do not have a significantly redder outer region, and 
their overall colors are bluer than Sbc's. 

Misclassifications: Two BSCs have red cores and blue outer 
regions (092_4957 and 294_2078). Both have tilted emission 
lines, indicative of rotation. Their rotation signature, as well as 
the overall structural and color gradient appearance suggest that 
they are more likely to be late-type galaxies rather than genuine 
BSCs. The remaining BSC (212_1030) seems to be an early- 
type galaxy with a red color similar to that of local early-type 
galaxies. 

Further evidence for the existence of red outer regions for 
Type-1 BSCs comes from the 2-dimensional bulge+disk SB 
profile fits. These provide the structural parameters and colors 
of photo-bulge and photo-disk components, which can be com- 
pared with the results obtained from the color gradient analy- 
sis. Since we are interested in obtaining well-constrained colors 
and structural parameters for each component, we use the "si- 
multaneous" fitting method, i.e., where we fit the /-band and 
V-band images of each BSC simultaneously, allowing only the 
total flux, B/T, and the background correction to be different 
in each band (Simard et al. 2001; Koo et al. 2001). 

To test the robustness of our results to the assumptions for 
the adopted SB profiles, we tried two kinds of SB models - 
the 'V 1 ' 4 + disk model" which uses a r 1 / 4 law profile for the 
photo-bulge and an exponential profile for the photo-disk, and 
the "double exponential model" which adopts exponential pro- 
files for both photo-bulge and photo-disk components. Note 
that the derived structural parameter values from these "simul- 
taneous" fits (Tables 2 and 3) can be slightly different from the 
values derived from the "separate" fits described in section 2, 
which is used for the sample selection (Table 1). 

Fig. shows the azimuthally averaged major axis surface 
brightness profile of BSCs (points) obtained from the aper- 
ture photometry, together with the best-fit model profile from 
GIM2D (solid line) which is composed of a photo-disk (dotted 
line) and a photo-bulge (dashed line). Fig. [7] shows that the 
surface brightness profiles cannot be fitted well with a single 
exponential profile - a significant second component is needed 
to fit their SB profile. In Tables 2 and 3, we list the GIM2D 
parameters from the "simultaneous" r 1,/4 +disk and double ex- 
ponential model fits, respectively. The important quantity in 
Tables 2 and 3 is the color of the more extended component, 

1. e., a component whose half-light radius is the larger of the two 
(fhi.maj = r e .fnaj for a r 1 / 4 component, and r hLmaj = 1 .68 r djnaj for 
an exponential component). Regardless of the assumed model 
SB profile, the extended components (typically, photo-bulge) 
of three Type-1 "good" BSCs indeed have a color redder than 
typical Sbcs (U - B ~ 0) or as red as the expected colors of 
passively-evolving E/SOs (U-B ;> 0.1), providing an indepen- 
dent confirmation for the existence of a red underlying stellar 



population in these BSCs. For the three remaining "possible" 
BSCs, the 2-D profile fits are less conclusive. Some redder 
components are not more extended. Nevertheless, even for such 
cases, the redder and bluer components have comparable ru,maj- 
We will use the term "underlying red stellar component" or "un- 
derlying red component" to describe the extended red compo- 
nent which seems to exist in Type-1 BSCs. 

In principle, the shape of the surface brightness profile of the 
underlying red stellar population can tell us about the likely lo- 
cal counterparts of these BSCs. If they are like massive E/SOs, 
the underlying stellar population will be distributed following 
the r 1 / 4 law profile. On the other hand, the underlying red 
component will show an exponential SB profile, if BSCs are 
more likely progenitors of less massive spheroids. Unfortu- 
nately, we find that the difference in reduced x 2 values from 
both fitting models is small (A\ 2 <J 0.01), meaning that the 

r 1 ' 4 +disk model and the double exponential model are almost 
equally good. Thus, the current data do not constrain the likely 
local counterpart of the BSCs in terms of the surface bright- 
ness profile. Any underlying old stellar population will stand 
out more clearly in the near-infrared (NIR; e.g., Hammer et al. 
2001); thus an analysis of high-resolution NIR images should 
offer more insights toward the characteristic SB profile and the 
mass content of the red component. 

4. DISCUSSION 

This section provides detailed discussions on our findings. 
In section 4.1, we discuss the robustness of the use of emis- 
sion line widths for mass estimates. Section 4.2 shows detailed 
comparison of our BSCs vs. other BSCs in the literature, and 
section 4.3 discusses BSCs in the Hubble Deep Field. A rough 
estimate of the number density of BSCs is given in section 4.4, 
and section 4.5 discusses BSCs and their interaction with close 
neighbors. 

4.1. Mass of the underlying stellar population 

In the previous section, we found that 5 or 6 Type-1 BSCs in 
our sample appear to have an underlying, extended, red stellar 
component (U-B > 0.0) and a bluer core from where most of 
the flux in emission lines may originate. When the star forma- 
tion is very localized, the line width from emission lines can 
be much narrower than a more global velocity dispersion. It is 
thus reasonable to ask how well the velocity dispersion mea- 
sured from the emission lines really represents the velocity dis- 
persion of the entire galaxy, especially for Type-1 BSCs. 

One possible way to resolve this issue is to directly com- 
pare a em with the observed absorption-line velocity dispersions 
(o~ a bs) of underlying stellar populations. Kobulnicky & Geb- 
hardt (2000) measure both o~ a t, s and c em for a sample of lo- 
cal galaxies of various morphological types, and find that a a b s 
and a an correlate well. However, their sample does not con- 
tain E/S0 galaxies, and thus may not be directly applicable to 
BSCs. For a sample of nearby elliptical galaxies, Phillips et al. 
(1986) measured er fm (core) from [Nil]. We identify 21 ellipti- 
cals in Phillips et al. (1986) which overlap with the elliptical 
galaxy sample of Faber et al. (1989) for which (7 fl /, s (core) is 
available. Here, we use the notation "(core)", since both <j em 
and a a b s are measured within r <J 0.1 r/,/. The comparison be- 
tween (7 fl /, s (core) and (7<, m (core) reveals that of these 21 mas- 
sive ellipticals, all have ^^(core) > 130 km sec -1 and none 
have <7 em (core)< 100 km sec" 1 . Some ellipticals in Phillips et 
al. (1986) do have er fm (core) < 80 km sec" 1 , but these galax- 
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ies are found to be rather less luminous systems (Mb < -19.5) 
and have [NII]A6584/// a ratios <; 0.2, resembling HII galaxies. 
Some BSCs have a em and Mb values very similar to the local 
HH-like ellipticals. A weakness in this comparison is that a are 
measured only near the core in Phillips et al. (1986), while a em 
of our BSCs are measured over a more extended region (~ r^if. 

For distant early-type galaxies, Koo et al. (2001) discuss 
a em vs dabs of luminous bulges at z ~ 0.8. They find 5 lumi- 
nous bulges for which a em is measurable, and find a reasonable 
agreement between a em and <r a t, s . None of these five objects 
have <r em < 100 km sec -1 . We find a similar result for the E/S0 
galaxy sample of Im et al. (2001), which contains at least 5 
objects with both a a b s and <j em . 

We would like to emphasize that none of the luminous bulges 
or E/SOs in either the nearby or distant samples have <j em < 100 
km sec" 1 . The lack of massive local ellipticals with a em < 80 
km sec" 1 suggests that Type-1 or Type-2 BSCs are more likely 
to be low-mass ellipticals with intense star formation rather 
than massive E/SOs. 

A more direct test would be to compare both a em and a a b s 
of the underlying old stellar population of BSCs. Unfortu- 
nately, our BSCs are in general too faint to provide contin- 
uum spectra with enough S/N for the measurement of a a bs- 
For this reason, we derive a limit on a„bs for only one object 
(103_6061) in our sample. Fig. || shows a continuum-divided, 
absorption-line spectrum of 103_6061, from which we obtain 
an upper limit a a b s < 80 km sec" 1 . The observed spectrum is 
fitted with a model spectrum using a maximum-penalized like- 
lihood method (Gebhardt et al. 2001; also see Saha & Williams 
1994; Merritt 1997). The model spectrum is comprised of high- 
resolution, stellar templates which are convolved according to 
a given velocity dispersion and instrument resolution and are 
mixed with proper weights to provide the best match to the 
galaxy spectrum. For more details on the absorption line fitting 
method, see Gebhardt et al. (2001). Some absorption lines in 
this particular object are contaminated by emission lines (dotted 
line in Fig. ||), and we artificially remove these lines when fit- 
ting the spectrum. In the figure, Balmer emission lines seem 
narrower than the stellar absorption lines, but this is simply 
because stellar lines are intrinsically broader than the nebular 
emission lines in the absence of velocity fields. The upper limit 
on cr a bs (80 km sec" 1 ) for 103_6061 is perfectly consistent with 
its a em value of 61 km sec" 1 , supporting the assumption that a em 
reflects the galaxy mass, or at least that <j em does not severely 
underestimate the kinematics of the galaxy. 

As an alternative check, we use the fundamental plane rela- 
tion of field, early-type galaxies at z <J 1 (Gebhardt et al. 2001) 
to estimate the velocity dispersion of the red, underlying com- 
ponents (hereafter, <jfp). This approach assumes that BSCs 
with underlying red, stellar populations, as defined in the previ- 
ous section, are early-type galaxies with a sprinkling of young 
stars at their center. The FP relation in Gebhardt et al. (2001) 
may be written as, 



logio(<T FP (kms 1 )) = 0.8 [logio(r h i M „ der (kpc))+032(SB MMnder )+9 

(2) 

where r„i has been replaced by rudder, the circular aperture 
half-light radius of the underlying, red component. SBid^nder 
is the rest-frame B-band surface brightness within ru, U nder of 
the same component corrected for the luminosity evolution de- 
rived in Gebhardt et al. (2001). Since there is uncertainty re- 



garding the plausible surface brightness profile for BSCs, we 
use structural parameters from both the simultaneous r 1,/4 +disk 
fit (Table 2) and the double exponential fit (Table 3), and use 
them to derive independent a F p values. The quantity ru^nder 
is derived as rounder = V^~ e b x r e,maj from Table 2 when the 
underlying red component is fitted with a r 1 / 4 law. When the 
underlying red component is fitted with an exponential profile, 
we use rhi^nder = 1 -68 ^/ cos(inc) rd t , na j from Table 2 or Table 3. 

In Table 4, we compare these a F p versus <r em . Note that 
the numbers within parentheses are the values derived from a 
double exponential profile. First, we find that a F p using the 
r 1 ' 4 +disk fit differ from a F p derived using the double exponen- 
tial fit. In many cases, the difference between two a F p values is 
small (±10 - 20%), but in some cases, a F p values from these 
different methods differ almost by a factor of two (172_5049 
and 103_6061). 

Second, the comparison between the a values shows that a em 
values are similar to or smaller than a F p by 20 ± 20% on av- 
erage. The derived ratio, — , is consistent with similar values 
derived from the comparison of velocity widths measured in 
HI, optical emission lines, and optical absorption lines for lo- 
cal galaxies (Kobulnicky & Gebhardt 2000; Telles & Terlevich 
1993; Rix et al. 1997). The amount of discrepancy is small 
enough to ensure that mass of BSCs in Table 1 should not be 
underestimated by more than a factor of ^2. In Fig. ||, we 
show ru unc ier vs. o F p of BSCs in Table 4. Like Fig. || we also 
plot the same relation for distant, red, QS-E/SOs (squares). The 
diamonds show rounder vs °~fp for the r 1 ' 4 + disk fit, and the 
stars indicate the same relation for the double exponential fit. 
Fig. H shows that BSCs still lie in the region which is populated 
by low-mass Es and dEs (Md yn < a f ew x 10 I0 M Q ) even though 
we use <r F p. A possible exception is 092_1339 (#1), which may 
have a dynamical mass of ~ 10 n M Q . However, even this ob- 
ject does not seem to be a typical massive E/S0 with > 10 1 1 M . 

In summary, the a em or "emission line widths" may not be a 
precise indicator of the kinematics of BSCs, but they seem to 
be useful as a rough estimator of the dynamical mass of BSCs 
(within a factor of ^2 uncertainty). All the different sources 
of evidence presented here (a em vs <r a b s of Es and a BSC, and 
a em vs a F p of BSCs) support our conclusion in section 3.1 that 
Type- 1 and Type-2 BSCs in our sample are likely to be either 
progenitors of dwarf galaxies, or of low- or at most moderate- 
mass E/SOs with Mdyn < a few x 1O 1O M . However, this con- 
clusion is based on relatively small number of objects, and more 
extensive studies are desired to establish the link between a em 
and Mdyn- 



4.2. Comparison with other studies 

Koo et al. (2001) study 68 high redshift (z ~ 0.8) luminous 
photo-bulges (Jbuige < 23.5) in the Groth Strip; five of their ob- 
jects overlap with ours. They find some "blue" photo-bulges, 
but these generally have relatively low luminosities. Since 
they will fade after the star-forming phase (e.g., Guzman et al. 
0621)996), such blue photo-bulges will not evolve into the more lu- 
minous bulges today. Their result is consistent with our findings 
that BSCs are relatively low mass objects. 

In order to check for a possible link between CGs and BSCs, 
we have also looked at values of B/T, R, and V -I for CGs in 
Phillips et al. (1997). Among 25 CGs at / < 22 with known 
z, we find only 2 that qualify as BSCs. The remaining 23 CGs 



When <7 em is measured over more extended region as for our BSCs, they are expected to be better tracers of the global kinematics. 
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turn out to be normal, red E/SOs (10), smooth galaxies with ex- 
ponential profiles (6), or galaxies with disturbed morphologies 
(7). 

Schade et al. (1999) studied elliptical galaxies in a subset 
of the CFRS and LDSS samples for which HST images are 
available. Eight ellipticals in the Schade et al. sample exist 
in the Groth Strip. Among these, 5 objects we classify as red 
QS-E/SOs, and one object is in our BSC sample (092_1339). 
The remaining two objects are not classified as E/SOs by us 
(084_1138 and 093_3251). A visual inspection of these two 
galaxies shows that they have red bulges surrounded by lumpy, 
blue low surface brightness regions, suggesting that they may 
instead be spiral galaxies, or red, early-type galaxies which are 
accreting gaseous blue satellites. All three possible blue ellip- 
ticals in Schade et al. (1999) mentioned above are either low 
mass, star-forming systems or galaxies with non-smooth mor- 
phology. On the other hand, one "possible" BSC in our sam- 
ple, 062_6465, is not classified as an elliptical in Schade et al. 
(1999). We suspect that this object is classified as SO in their 
work, since its bulge fraction is B/T ~ 0.5 and they classify 
objects with B/T ~ 0.5 as SOs. Also, note that Brinchmann et 
al. (1998) visually classify this object as an E/SO. 

Menanteau et al. (1999; 2000) studied spheroidal galaxies 
(E or SO) with HST and find that about 25% of these spheroidal 
galaxies have blue colors. More interestingly, they find that 
these galaxies have a blue localized region and a red outer re- 
gion similar to Type-1 BSCs. Thus, their spheroids seem con- 
sistent with ours. A subset of their sample comes from Abra- 
ham et al. (1999), and in order to understand similarities and 
differences, we have made a comparison of the Abraham et al. 
(1999) sample vs. BSCs, which is given in the next subsection. 

4.3. BSCs in the Hubble Deep Field 

The Hubble Deep Field (HDF) offers a much deeper HST 
WFPC2 image than our Groth Strip images. We have exam- 
ined the HDF images for galaxies that fit our selection criteria, 
and find seven QS-E/SOs with I < 22. Of these, only one (hd4- 
241.1) is a BSC. In contrast, Bouwens et al. (1998) identify 1 1 
E/SOs (/ < 22) in the HDF. Of these, Abraham et al. (1999) 
identify 4 with signs of recent star formation (see Menanteau 
et al. 2000 for additional work on this subject). The BSC hd4- 
241.1 (z = 0.321) is also found as a blue E/S0 by Abraham et al. 
(1999), but their 3 other blue E/SOs are not identified as BSCs 
or QS-E/SOs in our study, due to R values which slightly exceed 
0.1. 

Three of their 4 blue E/SOs, including hd4-241.1, have ap- 
parent magnitudes and sizes comparable to our BSCs, suggest- 
ing they are probably not massive E/SOs. The remaining object 
in their list (hd2-251.0 at z ~ 1) is somewhat exceptional, in 
the sense that it is a bright object (/ ~ 20) at high redshift for 
which we do not find a similar counterpart in our BSC sam- 
ple. This object is also quite red and has a blue core. How- 
ever, we note that this object has a very broad Mgll emission 
line indicating the existence of an AGN (Phillips et al. 1997), 
and is also detected as a radio (Fomalont et al. 1997) and X- 
ray source (Hornschemeier et al. 2000). Thus, the blue nu- 
clear color of this object probably originates from the AGN. On 
the other hand, the a measurement from absorption lines gives 
a ~ 195 km/sec (Gebhardt et al. 2001); thus the underlying old 
stellar population seems to be a massive system. We did not 
classify this object as an QS-E/S0 due to its slightly disturbed 

'The redshifts of companions are Zcomp = 1 -55 for 092_1339 (z = 0.90), Zcan 
Zcom,, = 0.88 (left) and Zcomp = 0.24 (right) for 212_1030 (z = 0.88). 



appearance (R = 0. 13) showing circular ripples of the type that 
can occur in merger remnants (e.g., Hernquist & Quinn 1987). 
In summary, our BSC criteria demand greater smoothness and 
symmetry than Menanteau et al. (1999, 2000) or Abraham et 
al. (1999), who find roughly twice as many blue spheroids in 
total. However, their blue E/SOs seem otherwise to have similar 
structural and photometric properties to our BSCs. Many even 
have blue localized regions and red envelopes like our Type-1 
BSCs. 

4.4. Number density of blue massive early-type galaxies 

We find that - 13 — 20% of QS-E/SOs with / < 22 are BSCs, 
depending on selection criteria. Among them, we estimate that 
none or perhaps only one (092_1339) out of 10 BSCs may be- 
come a massive early-type galaxy today with a ^> 150 km/sec 

or Md r n J> 10 n M Q . Thus, the fraction of truly massive blue 
E/SOs is found to be at most a few percent or less of the QS- 
E/S0 sample. This fraction might be larger since a star forming 
massive early-type galaxy shortly after merging may appear pe- 
culiar, thus falling outside our selection criteria. Studies of red, 
field E/SOs show that their number density at z < 1 is not much 
different from that of present-day E/SOs, providing an indepen- 
dent constraint on the fraction of blue merger remnants: the 
current uncertainty on the number density estimate for red, field 
E/SOs at z <^ 1 is about 30%, therefore the possible fraction of 
blue merger remnants is <; 30% of E/SOs assuming that merger 
events tend to create E/SOs than destroying/transforming them 
(e.g., Im et al. 2001,1999, 1996; Driver et al. 1998). 

4.5. Close neighbors 

Barton et al. (2000) studied close pairs of local galaxies 
(separation, Ar < 50 kpc), and find that star forming, bulge- 
like systems are possibly accreting gas from their close neigh- 
bors. By searching for close neighbors, we can test this sce- 
nario. The images of BSCs in Fig. ||have a physical scale of 50 
kpc by 50 kpc. Four BSCs have companion galaxies with simi- 
lar apparent magnitudes (092_1339, 103_6061, 172_5049, and 
212_1030). Redshifts of all these companions are available 9 , 
which show that the apparent companions of three of the BSCs 
are not physically associated with the BSCs. The only BSC 
with a physically associated companion within Ar < 25 kpc is 
212_1030, which is probably not a genuine BSC as we argued 
in the section 3.1. Thus, if the blue colors of BSCs are due to 
accreted gas, the entire companion was swallowed, leaving no 
discernible remnants. 

5. CONCLUSION 

From our DEEP GSS data containing 262 galaxies with 
/ < 22, we have identified 10 blue spheroid candidates at z <J 1, 
which we define as blue, bulge-dominated (B/T > 0.4), smooth 
galaxies morphologically similar to local E/S0 galaxies. We 
find that their internal velocity widths as measured from opti- 
cal emission lines are small, a em <J 80 km sec" 1 . Given their 
photometric, structural and kinematic properties, most BSCs 
resemble local H II region-like galaxies (e.g., Phillips et al. 
1986), and appear to be the progenitors of low-mass spheroidal 
galaxies today (Mdyn ^ 10 10 a fewM Q ). Two BSCs appear to 
be misclassified spiral galaxies, and one BSC seems to be a 
red early-type galaxy with blue close neighbors. We identify 

, p = 0.65 for 103_6061 (z = 0.36), z comp = 1.0290 for 172_5049 (z = 0.36), and 
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possible underlying red stellar components in 5 or 6 out of 10 
BSCs ("Type-1"). Velocity dispersions of the underlying red 
stellar component estimated from the FP relation or from a di- 
rect measurement of the absorption-line spectrum do not differ 
from a em values by more than a factor of 1.2 on average. This 
implies that even though the <j em values may underestimate the 
true kinematics of the galaxy, the difference is not enough to af- 
fect the conclusion that the underlying red components are not 
very massive (<; a few x 10 10 Mq). Overall, we find that mas- 
sive blue early-type galaxies with smooth, symmetric appear- 
ance are rare (a few percent of E/SOs with / < 22), suggesting 
either that the major merging events, which would produce star- 
forming massive, smooth, early-type galaxies, are rare at z < 1, 
or that the merger-triggered star formation activities subside by 
the time the merger products appear smooth and symmetric. 
Given our small sample size, we encourage future studies of 



more BSCs, to obtain a more conclusive result and to find any 
trend in the abundance and mass of BSCs as a function of red- 
shift. 

This paper is based on observations with the NASA/ESA 
Hubble Space Telescope, obtained at the Space Telescope Sci- 
ence Institute, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under NASA con- 
tract NAS5-26555. Funding for DEEP was provided by NSF 
grant AST-9529098. This work was also supported by the 
STScI grants GO-07895.02-96A, AR-.6402.01-95A, and AR- 
07532.01-96. We are grateful to Raja Guhathakurta, Garth 
Illingworth, Nicole Vogt, Vicki Sarajedini, and Andy Sheinis 
for their help in acquisition and reduction of spectra, and useful 
discussions. We also thank the anonymous referee for useful 
comments and suggestions. 
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FIG. 1. — Spectroscopic redshift vs. (V —I) color for 262 GSS galaxies with 16.5 < / < 22 (dots). QS-E/SOs are shown as open squares and triangles. Among 
them, BSCs are marked with thick triangles, and additional "possible" BSCs are marked with thin triangles (see text). The two dashed lines represent the upper and 
lower limit on colors of a passively evolving stellar population formed at Zf or = 11 (for more details, see text). BSCs are chosen as QS-E/SOs which do not meet this 
color criterion. 
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FIG. 2. — Postage stamp images of 10 BSCs (/-band). The first 6 images (from top left to right) are for "good" BSCs, and the remaining images are for "possible" 
BSCs. Parameters listed in Fig. □ are, from left to right, and from top to bottom, (1) ID# in our GSS catalog, (2) I to , magnitude, (3) redshift, (4) effective radius 
fv(arcsec), (5) velocity dispersion <r, (6) observed (V — /), (7) B-band absolute magnitude, (8) residual parameter R, (9) B/T. Boxes are scaled to be 50 kpc on its 
side. 
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FIG. 3. — Redshift vs. / magnitude of GSS galaxies. QS-E/SOs are marked with squares, "good" BSCs are marked with thick triangles, and "possible" BSCs 
are marked with thin triangles. The remaining galaxies are plotted as dots. Lines indicate trajectories of galaxies with two different absolute magnitudes at z = 0: 
M B (z = 0) = -19.37+5 log(A) (~ L* \ Marzke et al. 1998) and M B (z = 0) = -17.5 + 5 log(/i). The solid lines assume passive luminosity evolution equal to -1.7 X z in 
fi-magnitudes (Im et al. 2001; Gebhardt et al. 2001; Koo et al. 2001), while dashed lines do not take the evolutionary correction into account. For the ST-correction 
of the non-evolving models, we use the 13 Gyr-old SED of a model stellar population with 0. 1 Gyr burst SFR, z / or =11, and the Salpeter IMF (GISSEL96). 
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FIG. 4. — The Oil [3727] doublet for representative BSCs. Most spectra here come from LRIS data. For 092_1339, we also show ESI data. The comparison 
between the LRIS and ESI data supports the accuracy of the line width measurements from lower-resolution LRIS data. 
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FIG. 5. — a vs. Ru.cir for BSCs. The spaces occupied by typical local Hubble types are taken from Phillips et al. (1997). BSCs (triangles) are generally low-mass 
objects that lie in or near the space occupied by dEs). Red, QS-E/S0 galaxies with spectroscopic redshifts in Im et al. (2001) are plotted with squares showing a 
region where these QS-E/SOs with / < 22 are expected to be found. The number plotted next to each point is the object number that appears in the first column of 
Table 1. The dotted lines show loci of constant mass from Eq. (1). Object #8 is not plotted as we have no a measured for it. 
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FIG. 6. — Azimuthally averaged rest-frame U—B color profiles plotted along their major axes, and 2-dimensional U—B color maps of BSCs. For the 1-D color 
profile, the x-axis shows the distance from the center of the object, r, in 0. 1 arcsec. For the 2-D color map, the x and y axes are also in units of pixels. Type- 1 BSCs 
have a central region significantly bluer than the outer region, which in many cases is almost as red as typical field E/S0 galaxies. Note that the expected U—B colors 
of passively evolving stellar populations after 0.1 Gyr bursts range from 0.1 to 0.5 depending on the age and metallicity at the redshift of interest (see text). 
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FIG. 7. — Azimuthally averaged surface brightness profile of BSCs plotted along their major axis, as a function of distance from the center of the object. For 
each object, there are three panels, two for comparison with the r 1 / 4 bulge+ exponential disk profile, with distances indicated in units of r 25 (0. 1 arcsec) and r(0. 1 
arcsec), and a third for comparison with the double exponential profile, with distances given for r. A straight line is expected in the r 0,25 vs. SB plot for the r ' 25 
law profile, while the same is true in the r vs. SB plot for the exponential profile. The "X" marks show the measured SB. Models are drawn with the solid line 
(bulge+disk), the dashed line (bulge component), and dotted line (disk component), using outputs from the 2-dimensional GIM2D fit. Slight discrepancy between 
these data and the model exists; this is due to that (1) the aperture data are not deconvolved; (2) the aperture photometry is done with a fixed ellipticity while the 
bulge and disk components from the GIM2D fit have different ellipticities; and (3) faint neighbors affects the aperture photometry at the very outer regions. 
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FIG. 8. — The absorption-line spectrum of 103_6061 compared with the best-fit model spectrum. The spectrum is redshifted to the rest-frame and divided by the 
mean local continuum. The absorption-line velocity width (cr a ti S ) is measured using the maximum likelihood method as described in the text. We find cr a t s < 80 
km sec~' . The best-fit mix of stellar types is indicated in the figure. 
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FIG. 9. — Distribution of u FP vs. Rm, under f° r the underlying red components of Type- 1 BSCs. The spaces occupied by typical local Hubble types are taken from 
Phillips et al. (1997). Big diamonds indicate results based on a r 1 / 4 photo-bulge + exponential disk fit, while asterisks are for a double exponential fit. The line 
connecting diamonds and asterisks represents for each BSC how the position changes when the fitting method changes. The number plotted next to each point is 
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redshifts in Im et al. (2001) are plotted with squares as in Fig. 5. Even using parameters of the underlying stellar population, BSCs appear to be low-mass objects. 
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Table 1 

Basic Data of Blue Spheroid Candidates 
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0.87 
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±30 
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-0.04 


-0.03 


-0.04 




-0.1 








10 


273_4427* 


14:15:34.140 


52:05:53.96 
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21.56 


0.89 


-0.29 


-21.32 


0.40 


0.09 
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±30 
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NOTE. — (1) Object number as used in Fig. 5; (2) Source ID given by FFC-XXYY, where FF is the subfield, C is the WFPC2 chip number, and XX and YY are the chip coordinates in units of 
10 pixels; (3) Right Ascension (12000); (4) Declination (12000); (5) Spectroscopic redshift; (6) /-band total magnitude; (7) Observed V-I color derived from the total model-fit magnitudes; (8) 
Rest-frame U-B color, which is derived from V-/ and z using the A"-correction in Gebhardt et al. (2001); (9) Rest-frame B-band absolute magnitude; (10) Bulge to total light ratio measured in / 
band; (1 1) Residual parameter measured in /band; (12) Velocity dispersion a e „, measured from emission lines in kmsec -1 ; (13) Circular aperture half-light radius in pixels. One pixel corresponds 
to 0.1 arcsec; (14) Dynamical mass in units of 10 10 Mq. These mass estimates are uncertain by a factor of 2; (15) BSC Type as described in the text; Errors associated with each quantity are 
±1<j. 

'"Possible" BSCs. All others are "good" BSCs. 



Table 2 

Structural Parameters of BSCs from Simultaneous r 1 / 4 + DISK Fit 
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+0.2 


+0.1 
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-0.1 
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103_6061 


0.50 
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0.8 
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0.60 


0.26 


-0.46 


0.33 


12 
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+0.4 


+0.1 


+0.13 
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-0.3 


-0.1 


-0.14 
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4 


172_5049 


0.70 


1.6 


1.3* 


0.78 
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-0.28 


1.39 


0.27 


17 
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+0.03 


+0.1 


+0.1 


+0.05 


+0.26 


+0.05 


+0.37 


+0.02 
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-0.05 
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-0.05 
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5 


212_1030 


0.52 
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5.4 


1.79 


1.59 


0.19 


0.08 


0.02 


5 


1.051 






+0.04 


+0.2 


+0.3 


+0.20 


+0.20 


+0.12 


+0.12 


+0.02 
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-0.05 


-0.2 


-0.3 


-0.19 


-0.15 


-0.11 


-0.09 


-0.02 
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6 


294_2078 


0.46 


4.1 


3.6 


3.11 
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0.85 


-0.42 
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+0.11 


+1.2 


+0.2 


+0.77 


+0.22 


+0.32 


+0.14 


+0.08 
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-0.10 


-1.0 


-0.2 


-0.63 


-0.25 


-0.29 


-0.16 


-0.08 
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7 


062_6465 


0.51 


3.0* 


1.5 


3.59 


1.01 


0.88 


-0.20 
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42 


1.008 






+0.07 


+0.5 


+0.1 


+1.40 


+0.12 


+0.17 


+0.07 


+0.04 


+2 








-0.05 




-0.1 


-0.47 


-0.18 


-0.12 


-0.11 


-0.05 


-3 




8 


082_5252 


0.47 


10.6* 


2.1 


1.17 


0.67 


0.14 


-0.38 


0.07 


2 


1.102 






+0.01 


+0.2 


+0.1 


+0.04 


+0.03 


+0.05 


+0.03 


+0.01 


+i 








-0.01 


-0.2 


-0.1 


-0.05 


-0.04 


-0.05 


-0.03 


-0.02 


-l 




9 


153_6056 


0.33 


2.9 


1.5* 


0.40 


1.74 


-0.76 


0.24 


0.44 


50 


1.042 






+0.04 


+0.2 


+0.1 


+0.16 


+0.15 


+0.12 


+0.12 


+0.13 


+2 








-0.05 


-0.2 


-0.1 


-0.18 


-0.15 


-0.13 


-0.12 


-0.05 


-4 




10 


273_4427 


0.39 


3.3 


0.7 


1.01 


0.79 


-0.22 


-0.36 


0.41 


74 


1.083 






+0.04 


+0.7 


+0.1 


+0.19 


+0.09 


+0.11 


+0.06 


+0.08 


+i 








-0.03 


-0.9 


-0.1 


-0.14 


-0.11 


-0.09 


-0.07 


-0.24 


-l 





NOTE. — (1) Object number; (2) Object ID; (3) Bulge-to-total light ratio measured in the /-band from the simultaneous fit; (4) 
Major axis bulge effective radius in pixels. One pixel is equal to 0.1 arcsec; (5) Major axis disk scale length in pixels; (6) Observed 
V-I color of bulge; (7) Observed V-/ color of disk; (8) Rest-frame U-B color of bulge; (9) Rest-frame U-B color of disk; (10) 
Ellipticity of the bulge component; (11) Disk inclination angle in degrees. For a face-on system, this value is equal to 0; (12) Reduced 
X 2 value of the 2-D fit. 

a The circular aperture half-light radius, can be obtained as cir = ^/l-e^ X r e ma j or . 

b Major axis disk half light radius can be obtained as r M .d = 1.68 X r<f. ma jor, and the circular aperture disk half light radius is 
m,cir = Vco.s(Inc) X r hLi . 

*Components which represent the underlying, red, stellar population for Type-1 BSCs. 
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Table 3 

Structural Parameters of BSCs from Simultaneous Double Exponential Fit 



# 


ID 


B/T 


a b 

r dl,maj 


-i b 

^.maj 


(v-r, dl 


(V-I)J2 


{U-B) dhresl 


(U-B) d2 ,res, 




Inc L 


x 2 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


J 




U.JZ 


Z.V 


u.o 


1 ^"7 


l.uy 


u.uo 


-V.Z.J 




At. 


inn 






+0.05 


+0.4 


+0. 1 


+0.18 


+0.15 


+0. 1 1 


+0.09 


+0.06 


+4 








-0.04 


-0.4 


-0.1 


-0.17 


-0.15 


-0.10 


-0.09 


-0.06 


-7 




2 


HQ") ZLCK7 


U.J 1 


6 4 


2 4 


52 


1 08 


-0 43 


n 1 7 


4l 


^ ] 


i.UUJ 






+0.04 


+0.1 


+0. 1 


+0.12 


+0.15 


+0.12 


+0.19 


+0.03 


+2 








-0.05 


-0.4 


-0.2 


-0.12 


-0.16 


-0.11 


-0.18 


-0.03 






3 


103_6061 


0.71 


0.8 


2.7* 


0.75 


1.29 


-0.31 


0.26 


0.19 


26 


1.046 






+0.07 


+0.1 


+i.i 


+0.17 


+0.44 


+0.17 


+0.53 


+0.03 


+10 








-0.10 


-0.1 


-0.5 


-0.16 


-0.53 


-0.16 


-0.56 


-0.05 


-16 




4 


172_5049 


0.73 


1.5* 


0.2 


1.19 


0.79 


0.16 


-0.26 


0.21 


30 


1.033 






+0.02 


+0.1 


+0.1 


+0.08 


+0.17 


+0.09 


+0.17 


+0.03 


+« 








-0.04 


-0.1 


-0.1 


-0.07 


-0.14 


-0.08 


-0.14 


-0.03 


-9 




5 


212_1030 


0.78 


2.6 


0.1 


1.84 


1.74 


0.22 


0.16 


0.18 


80 


1.054 






+0.02 


+0.3 


+0.2 


+0.16 


+0.26 


+0.10 


+0.16 


+0.08 


+4 








-0.03 


-0.3 


-0.1 


-0.15 


-0.24 


-0.09 


-0.15 


-0.09 


-14 




6 


294_2078 


0.55 


2.2 


5.2 


1.97 


0.61 


0.30 


-0.52 


0.66 


45 


0.985 






+0.12 


+0.3 


+1.0 


+0.50 


+0.39 


+0.26 


+0.25 


+0.04 


+10 








-0.15 


-0.4 


-1.2 


-0.36 


-0.32 


-0.20 


-0.21 


-0.07 


-17 




7 


062_6465 


0.48 


2.0 


1.0* 


1.17 


1.93 


-0.11 


0.29 


0.19 


49 


1.016 






+0.21 


+0.7 


+0.1 


+0.66 


+0.80 


+0.35 


+0.34 


+0.08 


+n 








-0.27 


-0.5 


-0.2 


-0.62 


-0.88 


-0.38 


-0.48 


-0.10 


-10 




8 


082_5252 


0.34 


6.4* 


1.8 


0.81 


0.80 


-0.24 


-0.25 


0.14 


1 


1.101 






+0.02 


+0.1 


+0.1 


+0.06 


+0.03 


+0.06 


+0.04 


+0.02 


+3 








-0.01 


-0.2 


-0.1 


-0.06 


-0.04 


-0.06 


-0.04 


-0.02 


-1 




9 


153_6056 


0.48 


0.9 


2.6* 


1.00 


1.44 


-0.33 


0.01 


0.36 


55 


1.040 






+0.08 


+0.1 


+0.6 


+0.46 


+0.45 


+0.35 


+0.35 


+0.08 


+5 








-0.21 


-0.3 


-0.6 


-0.50 


-0.68 


-0.37 


-0.51 


-0.07 


-9 




10 


273_4427 


0.24 


2.4 


0.6 


0.90 


0.83 


-0.29 


-0.34 


0.39 


73 


1.083 






+0.04 


+0.1 


+0.1 


+0.28 


+0.09 


+0.17 


+0.06 


+0.08 


+3 








-0.05 


-0.2 


-0.1 


-0.24 


-0.09 


-0.15 


-0.05 


-0.10 


-2 





NOTE. — (1) Object number; (2) Object ID; (3) Bulge-to-total light ratio measured in the /-band from the simultaneous fit; (4) Major axis 
exponential scale length for exponential photo-bulge in pixels. One pixel is equal to 0.1 arcsec; (5) Major axis disk scale length in pixels; 
(6) Observed V -I color of exponential photo-bulge; (7) Observed V -I color of photo-disk; (8) Rest-frame U - B color of exponential 
photo-bulge; (9) Rest-frame U-B color of photo-disk; (10) Ellipticity of the exponential photo-bulge component; (11) Disk inclination 
angle in degrees. For a face-on system, this value is equal to 0; (12) Reduced \ 2 value of the 2-D fit. 



"The circular aperture half-light radius, r hl c „, can be obtained as r hl ar = y> 1 -e rfl(2 ) X (Vi(2), m aj. 

b Major axis disk half light radius can be obtained as r M .d = 1 .68 X r^ m aj . The circular aperture half-light radius is r M .dr = \/cos(Inc) X 
c ed2 = 1 -eoj(Inc). 

* Components which represent the underlying, red, stellar population for Type-1 BSCs. 



Table 4 

Comparison between Velocity Dispersion Estimates 



# 

(1) 


ID 

(2) 


<TFp(kms ') 
(3) 


a em 

(4) 


1 


092 1339 


136(92) ± 10 


0.57(0.85) ± 0.06 


3 


103 6061 


42(21) ±4 


1.44(2.87) ±0.73 


4 


172 5049 


111(61)±8 


0.53(0.96) ± 0.27 


7 


062 6465 


97(88) ±14 


0.29(0.33) ±0.31 


8 


082 5252 


72(50) ± 3 




9 


153 6056 


67(74) ±7 


1.18(1.06) ±0.46 



{^) a = 0.81(1.22) ±0.24 



NOTE. — (1) Object number; (2) Object ID; (3) Estimated a 
from fundamental plane relation for the underlying red com- 
ponent from the r 1 ^ 4 + disk fit and from the double exponen- 
tial fit (numbers in parentheses); (4) Ratio of a from emission 
lines i<J em ) to app 

a Mean and square-rooted variance of 



